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Abstract 

We show how to determine the ratio of the transverse velocity of 
a source to the velocity of emitted particles, using split-bin 
correlation functions. The technique is to measure S2 and S^, 



*On leave until 12 October 1993 at: Theory Division, CERN, CH-1211 Geneva 23, 
Switzerland; Internet: seibert@suryall.cern.ch. 
' Internet : haglin@hep . physics . mcgill. ca. 
* Internet : gale@hep . physics . mcgill. ca. 



1 



subtract the contributions from the single-particle distribution, and 
take the ratio as the bin size goes to zero. We demonstrate the 
technique for two cases: each source decays into two particles, and 
each source emits a large number of particles. 

PACS Indices: 25.70.Np, 12.38.Mh, 12.40.Ee, 24.60.Ky. 
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1 Introduction 

There has been considerable study of correlation functions in high energy 
and nuclear physics. However, with the exception of source size 
measurements using identical particle (Hanbury-Brown/Twiss) 
interferometry [|IJ, the connection between measurements of correlation 
functions and the physics of particle collisions is tenuous. For example, it is 
possible to measure the size of correlation sources (the number of pions 
produced) in high energy collisions |§, but no other source characteristics 
are currently identifiable. This is unfortunate, as the source size is 
anomalously large in ultra-relativistic nuclear collisions [Q, [3] (the sources 
decay into at least fifteen pions each ||), and the nature of these sources 
has not yet been determined. 

In this article, we propose to measure the velocities of correlation 
sources, the objects that make up the intermediate structure of collisions, 
by using the rapidity, y, and azimuthal angle, 0, correlation functions 



where measurements are made over a window from Yq to Yq + Ay. Here 




(1) 



and 




(2) 
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and 



cYq+ix pit rYo+ix r2n 

ay/x I dyJ M ab I dy 2 d^ h p {2) 



(2)Y ST jY o+d-l)x J JY +(i-l)x Jtt / j \ 

P h ^ xAY/4 ' { V 

where AY is an integer multiple of x. Split-bin correlation functions are 
useful for measuring correlations with maximum statistics (without re-use 
of data). S 2 is closely related to the standard two-particle correlation 
function: 

c (, Avn „ p (2) (0,fr/2) 

S 2 (Sy;AY) * p(2)(Q>Ay/2) • (5) 

For large AY, the correlation functions are almost independent of AY, but 
in practice it is impossible to make correlation measurements in particle 
collisions without using a large bin for reference, as the amount of produced 
matter fluctuates from event to event. 

The proposed technique is to measure the ratio 

1 ' S 2 (0;AY') - E 2 (AF) m 

Here 
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rYn + AY/2 pYnA-AY (7) 



and 
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Y<P(AY\ — -IS J JL 

A ^2\ 1 ^ 1 J r Y +AY riv rYo+AY r 2w ' 

/ dm / d^ h / dy 2 / #f p( Vl , p(y 2 , 0f) 

JY JO JY Jtt 

(8) 

are the values of S 2 and S 2 from the single-particle distributions, p(y) and 
p(y,4> lab ). In the case of a flat distribution [p(y, lab ) = const.], 



£ 2 = = 1. It is also possible to measure emission velocities using 
standard two-particle interferometry techniques; however, these 
measurements are somewhat more difficult than the technique proposed 
here, as they require particle identification. The technique proposed here 
has the advantage of being useable even in the absence of particle 
identification. 

We assume isotropic emission and show that, in the non-relativistic 
limit, R depends only on the velocity ratio, v s /v w , where v s is the source 
transverse velocity, and v n is the velocity of emitted particles in the source 
rest frame. In section |2|, we calculate R for the case of a resonance that 
decays isotropically into two particles. In section we calculate R for 
isotropic decay of a thermal source, finding that the shape is not very 
different from the first case. Finally, we summarize our results and discuss 
applications and future work in section £|. 



2 Two— particle decay 

Suppose that we have an event with a single resonance that decays into two 
particles, and a flat background of iV — 2 particles in a window of rapidity 
AY. The background contribution to the two-particle density is then 
constant: 

P ?\yiAi h ,y2A2 h ) = (N-2)(N-3)/(4n 2 AY 2 ) . (9) 



5 



Adding the contribution from combining one background particle with one 
particle from the resonance, we obtain the total uncorrelated signal, 

p^( yi ,^ b ,y 2 ,^ b ) = [iV(iV-l)-2]/(47r 2 Ar 2 ). (10) 

Given the resonance azimuthal angle, r , and the resonance rapidity, y r , 
we calculate the laboratory coordinates, y and lab , as functions of the 
center-of-momentum coordinates, and (f>. 

, / (\-vl) l l 2 v n cos0 \ 
V = tanh v . , + y r , 11 

y 1 + v s v tt sin & sin y 

i / + ^tt sin 6 1 sin \ > ,. , . . 

^ = tan' 1 " ■ n A + vriJ(-cos0) + r , (12) 

\(1 — t>s) ' fjr Sinti COS0 / 

where 

f x < 0, 
= " (13) 

[ 1 x > 0. 

[We use standard high-energy units, with c = 1. In deriving eq. flT2|), we 
assume that — 7r/2 < tan _1 x < tt/2.] We assume that y r and r are 
distributed randomly, with flat distributions. The resonance contribution to 
p^ is simplest in the resonance rest frame: 

mta a, a a\ sin0i sin0 2 5(cos0i + cos0 2 ) -<fo| - tt) 
Pr(^i,0i,^2,02) = ^ , (14) 

Z 7T 

where 5 is the Dirac 5-function. 

The resonance contribution to S 2 (5y; Ay) is 

rSy/2 rSy p2it p2tt 



S 2 , T (Sy;AY) 



Sy/2 JO 
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Sy*^AY>£\ yi , y2 ) 

(15) 



where 



p£ 2) (2/i>0i>2/2,02) 



80! 




d9 2 


dyi 




dy 2 



pi 2 \9 1 ,<f> 1 ,9 2 , ( f> 2 ) . 



(16) 



As we assume flat distributions, the contribution from all bins is the same, 
so we use only the bin from y — to 8y, and take Yq = The contribution 
at 5y = is 

" r8y/2 r2y r -8y/2 rZ8y/4 rSy/2 

/ dy r / dyi + dy r / 

/5y/4 jO J5y/2 J2y r -l 



4AY 



S 2 ,(0; AY) = lim ■ 



5y 2 [iV(iV-l) - 2] [(1-^ 2 )VV 

Ay 



2 [N (N — 1) — 2] (l-v*) 1 '^ 
Evaluating S 2l in the same manner, we obtain 

AY [ 2 d(j) r r dfa d(f) 2 5 (|0i - 02 1 - tt) 

q <f> (n. AY ) = Jo M ab =0 J 4>2 h =* 

2 ' rl ' J 2 7T 2 [iV(iV - 1) - 2] (1 - V 2)l/2 ^ 

Note that lab is now independent of 0, as in the limit — > the only- 
contribution is from cos0 = 0. We then perform the integral over r . 

Ay^d0s(|0 lab (0) - lab (0 + 7r)|) 



L , (17) 
.(18) 



(19) 



2tt 2 [N{N -1) -2]{1-v 2 s ) 1 / 2 v 7T 



where 



X < X < 7T , 

27T — X 7T < x < 2n 



g(x) = 
In the limit AY — > oo, 

<->2 — ^2 + <->2,r , ^2 2 "->2, 

Our ratio is then 



R 



1 r 2n 



2 , ^fl 1 



dab 



(0) - lab (0 + 7T)|) , 



(20) 



(21) 



(22) 



(23) 



Note that this result is independent of the size of the background, N\ In 
fact, if we allow the presence of more than one resonance, we find that S^r 
and S$ r are both proportional to the number of sources, n s , so R is also 
independent of the number of resonances. 

For the case of non-relativistic resonance decay, eq. fll2|) reduces to 

<T = tan- ( ^K)+ Bin0 \ + ^ ^ 
V cos © / 



as we can take COS6 1 = 0. It is clear from eqs. (PBQ and (f^j) that, in the 
non-relativistic limit, i? depends only on the ratio v s j v n . In Fig. 1, we 
show R for v n = c/2 and tv = c, along with the value in the non-relativistic 
limit. 

In both the relativistic and non-relativistic cases, -R(O) = 2; this is an 
identity, due to the different topologies of and y (0 = 2tc is equivalent to 
= 0, while y = Sy is not equivalent to y = 0). The behavior as — > 1 is 
analytically intractable, but from Fig. 1 it is clearly not too different from 
the non-relativistic case. The limit v s — > 1 can be taken analytically, 
yielding 



2 < 1 -f /2 l n 

The limit of large v s differs only by a factor 7, = (1 - v s 2 )i/2 in t he 
relativistic and non-relativistic cases. 



3 Thermal sources 



A typical example of a source that emits a large number of particles is a 
thermal source, such as a large droplet. Imagine that we have a single 
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source, as before, but that the source emits particles, where n n ^ 1. 
The two-particle density from this source is 

n n (n n — 1) sin^! sin0 2 



16 7T 2 



(26) 



The (^-function constraints are absent, as for large we can essentially 
ignore conservation of energy and momentum. 

Again moving to the resonance rest frame, and taking AY — > 00, we 
obtain 
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and 
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Taking the ratio as before, we find 
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(29) 



Note that R is independent of n w as well as being independent of N\ In 
fact, R is even independent of the number of sources, as before. 

We cannot evaluate the relativistic denominator for the thermal source, 
but we can take the non-relativistic limit as before. We obtain 



R = ^3 jTdfl sin fjd^ j\ 2 g{\^ h - 2 ab |) , 



(30) 



where 



,i a h 1 f(v s /v n ) + sin0sin<A rr . ,. , . . 

lab = tan" 1 v g/ ^ n )+ttH -cos0 ) + (j) r - 31 

\ sin cos / 



It is again clear that, in the non-relativistic limit, R depends only on v s / v n . 

In Fig. 2, we display R as in the case of two-particle decay. We find 
now that R(0) = 1 and R(oo) = 0, so the normalization is clearly different 
from that obtained for two-particle decay. Because the curves are different, 
there will be some dependence of the measured velocity on the source size. 
If the sources are known to emit many pions, as is the case in 
ultra-relativistic nuclear collisions, this dependence is weak. However, any 
velocity measurement using the proposed technique depends somewhat on 
knowledge of the source size. 

4 Conclusions 

We have demonstrated a new procedure for measuring the ratio of a source 
transverse velocity to the velocity of emitted particles, using the split-bin 
correlation functions S 2 and S$. We showed the technique for two simple 
cases: two-body decay, and unconstrained decay. In both cases, the 
technique gives a measurement that is independent of the number of 
sources, and the background size. In addition, we find that the 
measurement is independent of the source size for large sources. We have 
considered somewhat simple cases for pedagogical purposes, in order to 
illustrate our technique in a simple, model-independent framework. The 
specific issue of contamination by dynamical effects should be addressed 
separately. 

Similar techniques could be used to measure many other quantities, so 
the potential usefulness of this line of enquiry is enormous. For example, 
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comparison of energy and particle-number correlation functions could 
provide a relatively simple way to measure the energy distribution of 
minijets in ultra-relativistic nuclear collisions. These techniques are 
probably of the most use in high multiplicity collisions, where individual 
sources often cannot be resolved. We hope that our results will encourage 
other researchers to investigate new uses for correlation functions. 
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Figure captions 

1. R vs. tanh (v s ^ s /v n ) for two-particle decay. 

2. R vs. tanh (v s ^ s /v n ) for a thermal source. 
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Figure 1 
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Figure 2 



